Self-field effects upon the critical current density of flat superconducting strips 
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We develop a general theory to account self-consistently for self-field effects upon the average 
transport critical current density J c of a flat type-II superconducting strip in the mixed state when 
the bulk pinning is characterized by a field- dependent depinning critical current density J P (B), 
where B is the local magnetic flux density. We first consider the possibility of both bulk and edge- 
pinning contributions but conclude that bulk pinning dominates over geometrical edge-barrier effects 
in state-of-the-art YBCO films and prototype second-generation coated conductors. We apply our 
theory using the Kim model, J p k{B) = J p k(0)/(1 + \B\/Bq), as an example. We calculate J c (B a ) 
as a function of a perpendicular applied magnetic induction B a and show how J c (B a ) is related to 
J p k(B). We find that J c (B a ) is very nearly equal to J P K{B a ) when B a > B^, where B^ is the 
value of B a that makes the net flux density zero at the strip's edge. However, J c {B a ) is suppressed 
relative to J P K(B a ) at low fields when B a < B*, with the largest suppression occurring when B a /Bo 
is of order unity or larger. 

PACS numbers: 74.78.-w,74.25.Ha,74.25.0p 
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I. INTRODUCTION 



One of the most important physical properties char- 
acterizing a type-II superconductor is its transport crit- 
ical current I c . When the current / in a superconduct- 
ing strip exceeds I c , a voltage appears along the length; 
this voltage is due to the motion of vortices or antivor- 
tices across the strip. Applying a magnetic induction 
B a perpendicular to the current generally decreases I c , 
and another important physical property is the strength 
of the field dependence of I c (B a ) vs B a . For applica- 
tions it is of course desirable that this field dependence 
be as weak as possible, so that I c remains large in high 
fields. In low applied fields, the self-field generated by 
the current through the superconductor also suppresses 
the critical current. It is the purpose of this paper to 
present a method by which such self-field effects can be 
analyzed and calculated. 

The geometry of interest here is that of a long type- 
II superconducting strip of rectangular cross section 
sketched in Fig. 1. The strip is centered on the z axis 
and is of width 2a along the x direction and thickness 
2b along the y direction. For simplicity we assume that 
the thickness is somewhat larger than the London weak- 
field penetration depth A, such that (a) the sample di- 
mensions determine the vortex nucleation conditions and 
(b) we can use macroscopic thermodynamics, which in- 
volve the magnetic induction B, a local average of the 
microscopic magnetic induction over a few intervortex 
spacings, and the corresponding thermodynamic mag- 
netic field H = V B F(B), where F(B) is the Helmholtz 
free energy density in the mixed state. 2 When a magnetic 
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FIG. 1: Superconducting strip of rectangular cross section 
considered in this paper. Also shown are contours of constant 
vector potential A, which correspond to field lines of B, when 
the strip carries a uniform current density. 



induction B a is applied perpendicular to the z axis, B 
and H have only x and y components, which depend 
only on the coordinates x and y. We focus our atten- 
tion on dimensions of greatest interest to those involved 
in research and development on second-generation coated 
conductors [highly oriented YBCO (YBa2Cu307_a) films 
on textured substrates]: YBCO thicknesses in the range 
0.2-4 /xm and widths in the range 0.1-10 mm, such that 
b < a. 

Our paper is organized as follows. In Sec. II, we 
discuss the complications arising when both bulk and 
edge pinning are present, but we show that bulk pinning 
dominates over geometrical edge-barrier effects in state- 
of-the-art YBCO films and prototype second-generation 
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coated conductors. In Sec. Ill, we discuss the general 
problem of how to calculate the average critical cur- 
rent density J c (B a ) in an applied magnetic induction 
B a when the local depinning critical current density J p 
depends upon the local magnetic induction B(x,y) = 
xB x (x,y) + yB y (x,y). In Sec. Ill A, we give some basic 
equations, and in Sec. Ill B, we outline, but do not apply, 
a general method requiring a two-dimensional (2D) grid. 
In Sec. Ill C, we develop a general method for calculating 
J c (B a ) that uses only a ID grid but requires the integral 
of 1/ J p (B x , By) with respect to B x . This method ac- 
counts for the spatial variation of B{x, y) and J P {B) over 
the entire cross section of the film. In Sec. Ill D, we apply 
this method for perpendicular applied fields B a = yB a 
using the Kim model, 1 J pK (B) = J pK {0)/(l + \B\/B ), 
which assumes that J p depends only upon the magnitude 
of B(x,y). We calculate J c (B a ) vs B a and the field and 
current distributions across the strip. We briefly summa- 
rize our findings in Sec. IV. 



II. BULK VS EDGE PINNING 

According to critical-state theory in the mixed state 
of typc-II superconductors, 2 at the critical current of 
a superconducting strip the magnitude of the current 
density J(x,y) — J z (x,y)z = V X H(x,y) in a re- 
gion containing vortices is locally equal to the critical 
depinning current density J P {B). J P {B) characterizes 
the bulk pinning force per unit length f p on a vor- 
tex but is defined to have dimensions of current den- 
sity: f p — J p <po, where 4> = h/2e is the superconduct- 
ing flux quantum. On the other hand, the actual cur- 
rent density (averaged over a few intervortex spacings) 
is j(x,y) = j z (x,y)z = V X B(x,y)/(i , and it is this 
current density that enters into the Biot-Savart expres- 
sion from which the self-field contribution to B can be 
calculated throughout all space. 

Depending upon the magnitude of J p , the magnctic- 
flux-density and current-density distributions in a su- 
perconducting strip can be surprising complicated be- 
cause of the influence of geometrical edge barriers, Bean- 
Livingston surface barriers, or surface-barrier effects. 3-20 

For small J p , the self-field critical current of a super- 
conducting strip is determined entirely by the geomet- 
rical edge barrier. For a sample in the Meissner state 
subjected to an increasing current / in the z direction in 
the absence of a perpendicular applied field, the current 
density in the strip averaged over the thickness is 17,21 
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and the y component of the self-field in the plane y = 

. n 1-7 on 



is 17 - 21 



H y (x, 0) = 0, \x\ < a, 



27r|x|-\/a 



\x\ > a, 



(2) 
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to good approximation except within about b of the edges 
at x — ±a. Irreversible vortex penetration occurs and the 
critical current is reached when \H y (±a, 0)| first reaches 
the barrier-penetration field H s . In the absence of a 
Bean-Livingston barrier, H s = H c i, the lower critical 
field, but in the presence of a Bean-Livingston barrier, 
H s may be much larger, even approaching the bulk ther- 
modynamic critical field H c . 3 ' 22 ~ 25 However, the value 
of H s and the process by which irreversible penetration 
first occurs depend sensitively upon the local structure 
of the edge. 9 ' 11-14 ' 17,18 ' 26-31 The edge-barrier critical cur- 
rent in zero external magnetic field I s q can be estimated 
by evaluating Eq. (3) at x w a + 6/2, which leads to the 
result I s o = 2it\fahH s obtained in Ref. 17 for films with 
rounded edges. 32 Note that when I = I s q, the magnitude 
of the self-field in the middle of the top and bottom sur- 
faces is \H x (0,±b)\ = I s0 /2Tra = y/b/aH s , which is much 
smaller than H s when b <C a, the case of interest here. 

When the current / rises slightly above I s o, vortices 
and antivortices nucleate on opposite sides of the strip 
and then experience a Lorentz driving force per unit 
length of magnitude fd(x) = j z (x)(f> , which is sharply 
peaked at the edges and reaches a minimum at the cen- 
ter of the strip [see Eq. (1)]. Let us consider the case 
for which the bulk-pinning current density J p is indepen- 
dent of B (the Bean-London model 33,34 ). When J p is less 
than I s o/2irab, the driving forces exceed the bulk-pinning 
forces for all x, and the nucleated vortices and antivor- 
tices are driven towards each other until they annihilate. 
We therefore see that the self-field critical current I c of 
the strip is governed entirely by the edge barrier, i.e., 
Ic = IsOi independent of the bulk- pinning J p , when the 
ratio p = I p /Iso is less than 2/ir, where I p = AabJ p is 
defined to be the bulk-pinning critical current in the ab- 
sence of the edge barrier. 

For larger values of J p , i.e., for p > 2/ir, the criti- 
cal current I c of the strip (where a voltage first appears 
along the length) is larger than the current I s o when vor- 
tices are first nucleated. As the current I slightly exceeds 
I s o, the large screening currents near the sample edges 
drive nucleating vortices and antivortices toward the cen- 
ter of the strip, but they stop where j z (x) = J p . As 
shown in Ref. 19, further increases of the current produce 
metastable magnetic-flux distributions initially charac- 
terized by five zones. Near the edges are two vortex- 
free zones where the local current density j z (x) > J p . 
Along the centerline is another vortex-free zone where 
j z (x) < J p . On the right side of the central vortex- free 
zone is a region containing vortices with magnetic flux 
density B y (x) > 0, in which j z (x) = J p ; and on the left 
side of the central vortex-free zone is a region contain- 
ing antivortices with magnetic flux density B y (x) < in 
which j z (x) = J p . The critical current I c is reached when 
/ and j z {x) are just large enough to drive the vortex- 
filled and antivortex-filled regions together, shrinking the 
central vortex-free zone to zero width, such that vortex- 
antivortcx annihilation can occur along the centerline. 

In zero applied field when p is just above 2/V, the 
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state of the strip when I is barely above the critical cur- 
rent I c therefore can be described qualitatively in terms 
of four zones as follows: Vortices and antivortices nu- 
cleate at opposite edges of the strip, where the magni- 
tude of the local magnetic field is equal to H s . Once 
nucleated, vortices (antivortices) move rapidly to the left 
(right) through otherwise vortex-free zones near the strip 
edges where the supercurrent density j z {x) exceeds the 
bulk zero- field dcpinning critical current density J p . The 
vortices (antivortices) then enter a vortex-filled zone on 
the right (left) side of the centerlinc, where the magnetic 
flux density B y > (B y < 0) and the local critical cur- 
rent density j z {x) is just barely above the dcpinning crit- 
ical current density J p . The critical current I c is larger 
than 7 s o but also larger than I p because j z {x) > J p in 
the vortex-free zones. Analytic expressions from which 
I c can be calculated are given in Refs. 19 and 20 for the 
case that J p is independent of B. 

For increasing values of J p or p = I p /I s0 , the vortex- 
free zones at I = I c become narrower, and they shrink 
to zero width when the magnitude of the J p -generated 
self- field at the sample edges becomes equal to H s . If J p 
is independent of B and b <C a, the self-field at (x, y) = 
(±o, 0) can be shown from Eq. (9) to be 

H y {±a,0) = ±(J p 6/7r)pn(2o/6) + 1] (4) 

with correction terms inside the brackets of order (b/a) 2 . 
The value of J p that makes H y (a,0) = H s is J ps = 
TrH s /b{ln(2a/b) + 1], and the corresponding value of p 
is 

= Lps_ = 2 s Ja~Jb 

Ps I s0 In(2o/&) + l' (> 

where I ps = AabJ ps . 

For larger values of the depinning critical current den- 
sity, J p > J ps or p > p s , the critical current is given by 
I c = I p to excellent approximation because J z (x) = J p 
throughout the entire width of the strip except very near 
the edges. Any surface-current contribution at the edges 
arising from a surface-barrier-induced discontinuity in H 
can be shown to be smaller than I c by a factor of order 
b/a. 

Such complicated magnetic-field and current-density 
distributions arising from the geometrical edge barrier 
are important for strips fabricated of low-pinning su- 
perconductors, such as Bi-2212. On the other hand, 
for state-of-the-art YBCO films and prototype second- 
generation coated conductors, the depinning critical cur- 
rent density is so high that the parameter p is well above 
p s , and the excess current due to the edge barrier is neg- 
ligible. For example, Rupich ct al. 35 recently reported 
measurements of a self- field critical current of I c = 112 A 
at 77 K for a coated conductor composite of length 1.25 m 
and width 1 cm, in which the superconductor was metal- 
organic-derived YBCO of thickness 0.8 /xm coated on a 
deformation-textured NiW alloy substrate buffered with 
Y 2 3 /YSZ/Ce02. We can estimate I s0 as follows. The 



authors of Rcf. 36 found from magnetization measure- 
ments, using a 1 — (T/T c ) 2 temperature dependence for 
H C (T), V2H C (0) = 1.4 x 10 4 Oe, k = 57, and T c = 93.9 K 
for a YBCO single crystal, which yields H c = 3.24 kOe 
= 2.6 xlO 5 A/m at T = 77 K. Taking the nucleation 
field H s to be the same as H c i and using the Ginzburg- 
Landau expression H cl = (ff c / K \/2)(lnK + 0.5), we ob- 
tain H s = H cl = 180 Oc = 1.5 x 10 4 A/m and I s0 = 4.1 
A. Taking I p = I c = 112 A, we obtain p = I p /I s o = 27, 
such that p 3> p s , where p s — 3.1 [Eq. (5)]. This indi- 
cates that the vortex-filled region in which the current 
density is at the depinning critical current density fills 
practically the entire cross section at the critical current 
of coated conductors such as that reported in Ref. 35. 

III. BULK PINNING 

For the remainder of this paper, we limit our atten- 
tion to materials in which the depinning critical current 
density is so high that geometrical edge barriers have a 
negligible effect upon the critical current. We now focus 
on the questions of how to calculate the the average crit- 
ical current density J c = I c /Aab in an applied magnetic 
induction B a and how to determine the magnetic-field 
and current-density distributions at J c when the local 
depinning critical current density J p has a significant de- 
pendence upon the local magnetic induction B. 

A. Basic equations 

Consider a type-II superconducting strip with a ge- 
ometry as sketched in Fig. 1. If the current density 
associated with the magnetic induction B — V X A 
is j(x,y) — (l//xo)V X B(x,y), the vector potential 
A(x,y) = zA z (x,y) generated by j(x,y) can be calcu- 
lated from 

A zr (x,y,a,b) = 
Mo f a f b 

~ 4^ J J 3z(u,v)\n{(x - u) 2 + (y - v) 2 }dudv, (6) 

where the subscript r refers to the rectangular cross sec- 
tion. If j = zj z and the sheet current density K z = 2bj z 
are constant over the cross section of the strip, the re- 
sulting vector potential and the x and y components of 
the magnetic induction are 

A zr (x, y, a, b) = (n Q K z /2ir)a zr (x 1 y, a, b), (7) 
B xr (x, y, a, b) = (n Q K z /2ir)b xr (x, y, a, 6), (8) 
B yr (x, y, a, b) = (n Q K z /2ir)b yr (x, y, a, 6), (9) 

where expressions for the functions a zr , b xr , b yr , and their 
partial derivatives are given in Appendix A. Shown in 
Fig. 1 are contours of constant A zr (x,y), which corre- 
spond to B r field lines. 
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In the limit of vanishing thickness of the strip of width 
2a, simpler results but with singularities at the edges 
can be obtained by replacing j z in Eq. (6) by j z (x,y) — 
K z (x)S(y), such that the vector potential becomes 

A zs (x,y,a) = J"" K z (u)ln[(x-u) 2 +y 2 }du, (10) 

where the subscript s is a reminder that the results apply 
to a strip in the limit of zero thickness. If K z is uniform 
over the width of the strip, the resulting vector potential 
and the x and y components of the magnetic induction 
are 



A zs (x,y,a) 
B xs (x,y,a) 
B ys (x,y,a) 



(li K z /2ir)a zs (x,y,a), (11) 
{^K z /2^)b xs {x,y,a), (12) 
(poK z /2n)b ys (x,y,a), (13) 

where expressions for the functions a zs , b xs , and b ys are 
given in Appendix B. 

We are interested in calculating the critical current I c 
as a function of the applied field when the magnetic in- 
duction B a is applied perpendicular to the z axis. If 
one equates the current density j z to the dcpinning crit- 
ical current density J P (B) in Eq. (6), which determines 
the vector potential A describing the self-fields, one must 
self-consistently determine the spatial variation of the lo- 
cal magnetic induction B = B a + V x A. In general, 
Jp(B) is not constant over the cross section of the sam- 
ple, and the critical current density J c (B a ) = I c /Aab is 
therefore not equal to J p (B a ). 



B. General procedure with a 2D grid and J P (B) 

The following general numerical procedure, which is 
similar to that used by previous authors 37,38 for the disk 
geometry, could be used to account for self-field effects 
when the depinning critical current density J P (B) de- 
pends upon both the magnitude and direction of the lo- 
cal magnetic induction B. Divide the rectangular cross 
section into N — N x N y current elements of dimensions 
Ax x Ay, where Ax = 2a/N x and Ay — 2b/N y , centered 
at the coordinates (x n ,y n ), where n = 1,2, ...,N. If N 
is large, the bulk-pinning critical current density J p {B n ) 
in element n is very nearly constant. The local mag- 
netic induction in element n can be expressed as B n = 
B(x n ,y n ) = B a + B se i f (x n ,y n ), where B se i f {x n ,y n ) = 
xB se if, x (xn, y n ) + yB se if, y (x n , y n ) can be obtained with 
the help of Eqs. (8) and (9) by summing the contributions 
from all current elements: 

P i \ _ ST Mo^m, , Ax Ay 

£>seif,x\X,y) — / j o xr \x x m ,y y m , , ), 



ra—l 



N 



2 ' 2 
(14) 



B S eif,y(x,y) — y ] b yr (x x m ,y y m , , ) 



2 ' 2 
(15) 



with K m — J p (B m )Ay. This results in 2N equations, 
which for each given bulk-pinning critical current density 
function J P (B) can be solved by iterative methods to 
obtain the x and y components of the N vectors B n . The 
average critical current density can then be calculated as 
follows: 



1 N 

J c (B a ) = — J p (B n ). 



(16) 



n=l 



C. General procedure for strips with a ID grid and 

J p (B X , By) 

Suppose the strip is subjected to a uniform applied 
magnetic induction B a = xB ax + yB ay . We consider the 
case in which the underlying depinning critical current 
density J p (B x ,B y ) is a known function of both the x and 
y components of the local magnetic induction B. For ex- 
ample, J P (B X , B y ) could be determined experimentally 
by (a) measuring J c in applied fields B a large enough 
that the effects of the self- field B se i f are negligible and 
(b) fitting to a model that allows extrapolation to small 
B. In calculating the remanent magnetization of disks in 
self-fields, the authors of Refs. 37 and 38 found for the 
disk geometry that the magnitude of B changed so sig- 
nificantly over the disk thickness that the depinning crit- 
ical current density at the top and bottom surfaces was 
significantly suppressed relative to that midway between 
the two surfaces. This suppression of the depinning cur- 
rent density near the top and bottom surfaces must also 
be taken into account to achieve self-consistent solutions 
accounting for self-field effects upon the average critical 
current density J c (B a ) vs B a of long strips. 

There are three properties of thin superconducting 
strips with k<a (see Fig. 1) that enable us to accomplish 
this using only a ID grid by applying an idea from Refs. 
38 and 39. We simplify the general procedure outlined in 
Sec. Ill B by taking N y = 1, such that N = N x , and by 
considering current elements of dimensions Ax = 2a /N 
and Ay = 2b, centered at the coordinates (x n ,0), where 
n = 1,2, N. For simplicity, we use the approximation 
that B w fM)H, such that J w j = (l//x )V X B. 

The first important property of thin strips we use is 
that the x and y components of B se if(x,y) can be ob- 
tained to good accuracy by writing 

B se if,x(x, y)=^2 ^ m b xr {x - x m , y, ^-,b), (17) 



m=l 
N 



Ax 



B S elf,y{x, V) = Y1 /i " 2 ^ m *V( X ~ V, ~if> b )> ( 18 ) 

m—1 

where K n = 2bj z (x n ) is the average sheet current density 
in the element at x = x n . The net magnetic induction is 
B = B a + B self , where 

(19) 
(20) 



B x (x,y) = B ax + B ae if tX (x,y), 
B y (x,y) = B ay + B se if , y (x,y). 
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As can be seen from Fig. 1, B se if tX (x,—b) > and 
B se if,x{x,b) = -B aelftX (x,-b). 

The second property of thin films we use is that, al- 
though the current density j(x,y) = j z (x,y)z is given 
in general by j z = (l/no)(dB y /dx — dB x /dy), it is a 
good approximation for most of the strip width to write 
j z = —(l/fio)dB x /dy when 6 <C a. The accuracy of this 
approximation can be estimated for a uniform current 
density in a strip of rectangular cross section by calcu- 
lating the ratio Rj(x,y) = (dB y /dx)/(—dB x /dy) using 
Eqs. (A6) and (A7). Although Rj ~ 1 at the edges of the 
strip, we find that i?j(0,0) w (2/n)(b/a) at the center of 
the strip when b/a < 1. When b/a = 0.01, Rj(0,0) = 
0.0064 and Rj(x,0) < 0.1 for a - \x\ > 3.56. When 
b/a = 0.001, Rj(0,0) = 0.00064 and Rj(x,0) < 0.1 for 
a — \x\ > 3.56. 

The third property of thin films we use is that 
\dB y /dy\ <C \dB x /dy\, such that although B x (x,y) 
changes considerably over the film thickness, B y (x,y) is 
very nearly a constant, equal to B y (x,0). The accuracy 
of this approximation can be estimated for a strip car- 
rying a uniform current density by calculating the ratio 
R B (x) = \[B v (x,b) - B y (x,0)}/B y (x,0)} using Eq. (A3). 
Although Rb{x) is largest near the edges of the strip, 
R B {x) < 0.01 for a - \x\ > 3.56 when b/a < 0.01. 

The second property allows us to obtain the aver- 
age of j z (x,y) over the sample thickness from j z (x) = 
B se if yX {x,—b)/n,Qb, and the third property allows us to 
set B y (x, y) = B y (x, 0) in the critical state force-balance 
expression 

jz(x,y) = J p (B x (x,y),B y (x,y)), (21) 

such that in Eqs. (17) and (18) we may write 

K n = 2bj z (x n ) = 2B selftX (x n , -b)/no, (22) 

where B se if^ x (x n ,—b) is obtained either numerically or 
analytically from (see Refs. 38 and 39) 



B ax +B selfiX (x n -b) 



du 



c ( Xn -b) J P {u,B y (x n ,0)) 



= 2/x 6. (23) 



One may now use the above equations to find self- 
consistent solutions for B and j. The 2N unknowns, 
B x {x n ,—b) and B y (x n ,0), can be obtained numerically 
from 2N equations obtained from Eqs. (18-20), (22), and 
(23). Once these solutions have been found, we may cal- 
culate the average critical current density from 

1 N 1 N 

Jc(B a ) = -ft^jzM = Jj- -Y, B selfAXn,-b)(24) 

™ 71=1 



n=l 



D. 



Calculations in a perpendicular field with a ID 
grid and a Kim-model J p (B(x,y)) 



We have applied the numerical procedure described in 
Sec. Ill C to calculate J c {B a ) in a perpendicular field 



B a = yB a for three ficticious YBCO samples (a, b, and 
c) with dimensions as in Fig. 1 but with a = 120 /im and 
6 = (a) 0.1 fim, (b) 0.5 /zm, and (c) 1.5 /im. Suenaga et 
al. 40 reported that the depinning critical current densities 
inferred from 77 K magnetization measurements on three 
YBCO disks with these thicknesses were well fit by the 
Kim model, 1 



J p k(B) = J pK (0)/(l + \B\/B ), 



(25) 



except at low fields where self-field effects are responsible 
for deviations. For J p (B x ,B y ) in Eq. (23) we therefore 
took the local depinning critical current density to be 

J p (B x ,B y ) = J pK (B), where B = ^bJTb^. The pa- 
rameters were (a) J p k(0) = 4.32 xlO 6 A/cm 2 and B 
= 17.6 mT, (b) J pK (0) = 3.22 xlO 6 A/cm 2 and B = 
18.8 mT, and (c) J pK {0) = 2.55 xlO 6 A/cm 2 and B 
= 22.6 mT, respectively. 40 With Eq. (25) giving the ex- 
plicit form for J p , it was possible to evaluate Eq. (23) 
analytically with the result 



2B B x (x n , -6) + B x (x n , -b)^B 2 x (x n , -6) + B 2 (x„, 0) 

+ B 2 y {x n , 0) sinh" 1 [B x (x n ,-b)/\B y (x n ,0)\] 

= 2 i i Q J pK {Q)B Q b. (26) 

Using Mathematical 1 we solved Eqs. (18-20), (22), 
and (26) numerically with N — 101 to obtain B x {x n , —6) 
and B y (x n ,0) for samples a, b, and c. This allowed us 
to calculate B x {x,y) and B y (x,y) from Eqs. (17), (18), 
(22), and 

B x (x,y) = B selfiX (x,y), (27) 
B y (x,y) = B a + B selLy (x,y). (28) 

We also calculated B*, the value of B a at which 
B y (— a,0) = 0, for the three samples considered in this 
paper: B*/B a = (a) 0.651, (b) 1.320, and (c) 1.834 or 
B* = (a) 11.5 mT, (b) 24.8 mT, and (c) 41.4 mT. Shown 
in Fig. 2 are values of J c {B a ) vs B a calculated from Eq. 
(24). When B a > B*, the calculated values of J c (B a ) 
agreed with J p k{Bo) within (a) 0.6%, (b) 1.3%, and (c) 
2.1% for the three samples, where the values of B* a are 
shown by the solid symbols. Note that the suppression 
of J c (0) is greatest for the thickest film. In general, self- 
field effects become more important as the strip thickness 
increases. 

Shown in Fig. 3 are normalized plots of B y (x,0) vs 
x/a for the thickest of the three samples for applied fields 
B a = (a) 0, (b) B*/2, (c) B*, and (d) 3B*/2. Note that 
B y = at the left edge of the sample when B a = B*. 
Shown in Fig. 4 are corresponding normalized plots of 



j z (x) = B x (x,-b)/fj, b 



(29) 



vs x/a. Note that j z {x) is largest at those values of x 
where B y (x,0) = in Fig. 3. Note also that when B a > 
B* a , j z vs x becomes flatter; this is another indication 
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FIG. 2: The normalized average critical current density 
J c {B a )/J pK (0) vs B a /Bo calculated from Eq. (24) for the 
three samples a, b, and c (solid), compared with the Kim 
model, J p i((B a )/ Jpn(O), Eq. (25) (dashed). The values of 
Ba/Bo for the three samples are identified by the (a) solid 
circle, (b) solid triangle, and (c) solid square. 

that self- field effects have less of an effect upon J c (B a ) 
as B a increases. 

A rough (within about a factor of 1.2) estimate of 
B*, perhaps useful in the analysis of experiments, is 
lioHy(a,0), where J p in Eq. (4) is replaced by the exper- 
imental value of J c (0). Using the solid curves in Fig. 2 to 
obtain such an estimate yields (a) 12.7 mT, (b) 30.1 mT, 
and (c) 49.6 mT for the three samples discussed above, 
for which J c (0)/J pK (0) = (a) 0.838, (b) 0.652, and (c) 
0.534, or J c (0) = (a) 3.62 xlO 6 A/cm 2 , (b) 2.10 xlO 6 
A/cm 2 , and (c) 1.36 xlO 6 A/cm 2 , respectively. 

Several features of the above behavior of J c vs B a using 
the Kim model 1 are close to those found analytically for 
J c (B a ) when the applied field is parallel to an infinite 
slab, as discussed in Appendix C. 



IV. SUMMARY 

In Sec. II, we briefly reviewed the behavior of the av- 
erage critical current density J c (B a ) of superconducting 
strips when both bulk and edge pinning play a role. The 
main conclusion of this section is that edge pinning has a 
negligible effect upon J c in state-of-the-art YBCO films 
and prototype second-generation coated conductors. 

In Sec. Ill, we considered the problem of how to ac- 
count self-consistently for the dependence of the local 
depinning critical current density upon the local mag- 
netic induction when calculating J c (B a ). In Sec. Ill A, 
we presented some basic equations, and in Sec. Ill B, we 




-3-2-10 1 2 3 

x / a 



FIG. 3: Normalized magnetic flux density B y (x,0) vs x/a for 
the thickest of the three samples, calculated from Eq. (28), for 
applied fields B a = (a) 0, (b) B*/2, (c) B*, and (d) 3B*/2 
(applied fields shown as horizontal dashed lines) , where B* = 
41.4 mT and B*/(fi bJ pK (0)/n) = 2.71. 

set up, but did not apply, a procedure that, using a 2D 
grid and starting from an underlying local depinning crit- 
ical current density J p {B x , B y ), would yield the average 
critical current density J c {B a ) = I c /Aab vs an applied 
magnetic induction B a for a sample of rectangular cross 
section of arbitrary width 2o and thickness 2b. 

In Sec. Ill C, we developed a general method that, 
using a ID grid and starting from an underlying local 
depinning critical current density J p (B x ,B y ), yields the 
average critical current density J c (B a ) = I c /Aab as a 
function of the magnetic induction B a applied to a flat 
sample whose width 2a is much larger than its thickness 
2b, a case of considerable interest to those involved in 
research on prototype second-generation YBCO coated 
conductors. This method should work well even when J p 
is a function of both the magnitude B and the direction 
B of the local magnetic induction B = BB. 

In Sec. Ill D, we applied the above method to calculate 
J c (B a ) in a perpendicular field B a using the Kim model 1 
[Eq. (25)], in which the local depinning critical current 
density J p is a function of the magnitude B but not the 
direction B of B = BB. The functional form of the 
B dependence of the Kim model's J P (B) resulted in the 
particular form of Eq. (26), which we used in determining 
B x (x, —b) and j z (x). However, our general method is not 
limited to the Kim model but could be used with any 
other realistic model (whether analytic or numerical) for 
J P (B). Instead of using Eq. (26), one could use Eq. (23) 
to obtain a new algorithm connecting B x (x n ,—b) and 
B y (x n ,0) from which J c (B a ), B y (x,0), and j z (x) could 
be determined. 
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constant) can be expressed as in Eq. (7), where 



a zr (x 1 y,a,b) = 



y + b. 



{(x — a) 2 [arctan( y ' - ) — arctan( 



x — a 



y-b, 



x — a 

,y + b, 



+ (x + a) [arctan( J ) — arctanf- 

x + a x + a 

,x — a 
y-b' 
x + a 



x + a, 



+ (y — b) [arctan( -) — arctan( 



(y + 6) 2 [arctan( 



y-b' 
x — a. 



arctan(- 



)] 



y + b y + b 

+ (x -a)(y + b) \n[(x - a) 2 + (y + b) 2 } 
-{x-a){y- 6) \n[(x - a) 2 + (y - b) 2 } 
+ (x + a)(y-b) ln[(.x + a) 2 + (y - b) 2 ] 
- (x + a)(y + b) ]n[(x + a) 2 + (y + 6) 2 ]}/46 



(Al) 



FIG. 4: Local current density averaged over the sample thick- 
ness jz(x)/J P K(0) vs x/a for the thickest sample, calculated 
from Eq. (29) for applied fields B a = (a) 0, (b) B*/2, (c) B*, 
and (d) 3B*/2, where B* = 41.4 mT. 



The x and y components of the magnetic induction de- 
rived from Eq. (7) and B — V X A are given in Eqs. (8) 
and (9), where 



In this paper we have considered only static current 
and field distributions at the critical current. For ex- 
tensive treatments of the dynamics of current and field 
penetration into strips, see Refs. 39,42-44. 



Acknowledgments 
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APPENDIX A: A z , B x , AND B y FOR A STRIP 
CARRYING UNIFORM SHEET CURRENT 
DENSITY 



b yr {x,y,a,b) = 

x-a y 
-|arctan(- 



2b 

x + a 



-) — arctan(-^ 



+ 



26 



[arctan( 



x — a 
x + a' 



x — a 



arctan( 



y-b , 

x + a' 



y-6 (x-a) 2 + (y-b) 2 ^ 
+ T ln[ " 



(x + a) 2 
(x + a) 2 



(y-b) 21 
(y + b) 2 , 



(x - a) 2 + (y + b)^ 



(A3) 



When a strip of rectangular cross section (width 2a 
and thickness 26) carries a uniform sheet current density 
K z = 2bj z as described in Sec. I, the vector potential 
derived from Eq. (6) (neglecting an unimportant additive 



Partial derivatives of the components of B, such as 
dB xr /dy — (noK z /2ir)db xr /dy, are given by the follow- 
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ing equations, 

db xr (x,y,a,b)/dx = 

(x + ay + (y + o) (X — ay + (y — by 
db yr (x, y, a, b)/dy = -db xr (x, y, a, b)/dx, (A5) 
db xr (x,y,a,b)/dy = 

[arctan( X — y ) + arctan(- — y) 
y-b y + b 

— arctan( — -j-^) — arctan( -)]/2b, (A6) 

db yr (x, y, a, b) / dx = 

[arctan( — — -) + arctan( - — - ) 

X CL X | CL 

- arctan(^-±-^) - arctan(^-^)]/2o. (A7) 



netic field is applied perpendicular to a flat strip, for 
comparison we present here analytic results for the av- 
erage critical current density J c (B a ) in the z direction 
when a magnetic induction B a is applied in the y direc- 
tion parallel to the surface of an infinite slab of thickness 
2a. We consider two models for the depinning critical 
current density J P {B) and assume that B = HqH, such 
that at the critical current we have for each model 



J z (x) = Jp{B y ) = 



1 dB y (x) 
Hq dx 



(CI) 



Associated with this current is the self- field B s , such that 
B y (±a) = B a ± B s . Integrating Eq. (CI) from one side 
of the slab to the other, we find that the average critical 
current density is 



Jr. 



Bs_ 



(C2) 



APPENDIX B: A z , B x , AND B y FOR A THIN 
STRIP CARRYING UNIFORM SHEET 
CURRENT DENSITY 

When a thin strip of width 2a and negligible thickness 
carries a uniform sheet current density K z as described in 
Sec. I, the vector potential derived from Eq. (10) (neglect- 
ing an unimportant additive constant) can be expressed 
as in Eq. (11), where 



a zs (x,y,a) — y[arctan( ) — arctan( X a )] 

V V 



+ ln[(x af +y*]- X -±^ \n[{x + af + y\ 



(Bl) 



The x and y components of the magnetic induction de- 
rived from Eq. (11) and B — V X A are given in Eqs. 
(12) and (13), where 



where the dependence of B s upon B a can be determined 
by solving Eq. (CI) for B y (x). For both models of J P (B), 
we find that J c (B a ) has an inflection point at B a = B* 
and is equal to J p (B a ) for \B a \ > B*. 



1. J v {B)=a/\B\ 

For the model J P (B) = a/\B\, we find that B* a = 
^/fio aa an d 



wi> = f- ( m^- m - B: - (C3) 



2. J P (B) = J p0 /(1 + \B\/B ) 



(C4) 



b xs (x, y, a) = arctan(^ — -) — arctan ^ — - )], (B2) 

y v 



For the Kim model 1 , J p (B) = J p0 /(1 + \B\/B ), we 
find that 



and 



h ( \ 1 ^ < ( x + a ) 2 + y\ 

b ys (x,y,a) = -H {x _ a)2+y2 ] 



(B3) 



and 



B* a = (B /2)Ul + AfioJpoa/Bo - 1) (C5) 



APPENDIX C: J c (B a ) FOR B a PARALLEL TO A 
SLAB 

Although our main interest in this paper is in the av- 
erage critical current density J c (B a ) vs B a when a mag- 



J c (B a ) VI + 2b* + 2b* 2 -b 2 -l 



pO 



b*(l + b*) 
-,b>b*, 



1 + 6' 

where b = \B a \/B and b* = B*/B . 



,0<6<&*(C6) 
(C7) 
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